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AUTHORS’  SUMMARY 


A  solution  has  been  found  for  the  problem  of  flow  about  a  profile  by  an  ideal 
incompressible  liquid  in  a  tunnel,  the  upper  and  lower  walls  of  which  may  have  varying 
porosity,  depending  on  whether  gas  flows  in  or  out  of  the  tunnel.  Calculations  have  been 
made  of  the  distributions  of  pressure  on  the  walls  of  the  tunnel  and  on  the  profile  located  at 


an  angle  of  attack  on  the  axis  of  symmetry  of  the  tunnel. 
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INTRODUCTION 

1  It  is  assumed  that  in  the  woiidng  section  of  the  tunnel  the  flow  is  incompressible  and 
inotational  or  potential,  and  that  the  upper  and  lower  perforated  walls  represent  a  tunnel  of 
infinite  length.  The  porosity  coefficients  of  the  upper  and  lower  wall  are  f+  and  f^ 
respectively,  outside  the  tunnel  is  surrounded  by  a  cavity,  in  which  the  pressure  is  constant 
(Fig  1).  The  profile  is  situated  on  the  mean  iii%  of  the  tunnel,  the  angle  of  inclination  of  the 
contour  of  the  profile  to  the  axis  of  the  abscissa  h±(x)  are  considered  to  be  small:  (h+fx)  - 
fOT  the  upper  part  of  the  profile,  and  h-(x)  -  for  the  lower  part). 

The  region  D  in  which  a  solution  of  the  problem  is  sought,  represents  the  interior 
of  a  two-dimensional  tunnel  having  a  width  H ,  with  a  slit  along  the  interval  [0, 1],  where  a 
slender  lifting  profile  is  located. 

In  a  linear  formulation,  by  transferring  the  bount^ary  conditions  to  the  axis  x  ,  wc 

obtain 

V(x,  +  0)  =  h+(x),  V(x,-0)  =  h-(x)  ,  (1) 

where  V  is  the  perturbed  vertical  component  of  velocity,  referred  to  the  velocity  of  the 
incident  flow.  The  conditions  on  the  walls 

y  =  y  .•x€(-oo,x),  V  =  RjU;  x6(Xj,oo),  V  = -R^U; 

JJ 

y  =  -y :  x6(-«>,X2),  V  =  R^U;  xefx^;***),  V  =  R“U  , 

where  U  is  the  perturbed  horizontal  component  of  the  component  of  velocity,  referred  to 
the  velocity  of  the  incident  flow,  xi  and  X2  are  the  coordinates  of  the  points  at  which  the 
porosity  is  equal  to  0 .  R^ ,  R^ ,  R” ,  R^  are  the  parameters  of  porosity  dependent  upon 
the  coefficient  of  porosity  f*"  and  f~ ,  and  also  on  whether  the  gas  flows  into  or  out  of  the 
tunnel.  Fig  2  gives  the  results  of  the  investigation  of  the  porosity  of  porous  panels  in  Ref  4. 

The  solution  of  the  problem  posed  is  reduced  to  finding  the  analytical  function 
<|»(z)  =  U  -  iV  of  the  complex  variable  z  =  x  +  iy  in  the  region  D  which  satisfies  the 
boundary  conditions  (1)  and  (2). 

2  Let  us  break  down  the  procedure  for  obtaining  a  solution  into  several  stages: 

(a)  Rnding  in  the  region  D  an  analytical  function  L(2)  =  LR-iLi  which  satisfies  the 
boundary  condition  (2)  solely  on  the  lines  y^  =  ±  ~  while  being  continuous  across  the 
sectitH)  xe  [0, 1] ,  y  =  0 . 
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Let  us  transfer  the  strip  in  the  plane  2,  bounded  by  the  lines  y±  =  ±  ^  .  by  the 
confOTinal  transformation  w  =  i  exp  z) ,  into  the  upper  half-plane  w  =  ^  +  ill  (Fig  3a). 
For  L(w(z))  the  conditions  (2)  on  the  true  axis  in  the  plane  w  assume  the  form 


^e(-«>o.^i(xi)). 

Li  =  -I^Lr; 

(3.1) 

(§i(xi),  0), 

Li  =-R;Lr; 

(3.2) 

^€(0,42(X2)), 

Li  =  RjLr; 

(3.3) 

(^2(X2).  *»)» 

Li  =  R^  Lr  . 

(3.4) 

Let  us  write  the  analytical  function  L(w)  in  a  general  form; 


I 
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L(w)  =  (w  -  (y/p  a  ^ 


where  a  is  a  cotuplex  constant,  P(w)  and  Q(w)  arc  arbitrary  polynomials  with  real 
cocfficicr  s  which  arc  independent  of  conditions  (3.1  to  3.4),  and  Yi,  Y2.  Y3  arc  real 


constants. 


Let  us  determine 


“R  “I  I  2 


fiom  the  conditirm  (3.4) 


a  =  exp(-i  arctg  R^)  ;  for  ^6  (0,  ^) ,  i)  =  0 


L(w)  =  Lj^-iLj  =  U  ’ 

from  the  condition  (3.3) 

Yj  =  i  (arctg arctg  R“) +nj ,  where  nj  =  1,0;  for  ^6(^j,0)  , 


TJ  =  0 


L(w)  =  Lj^-iLj  =  U  -^jrM^I^^U-^r‘fli^cxp(ijryj+i7ry2  + 

Q(^ 


m 


+  i  arctg  R^)  . 


From  the  condition  (3.2) 


Yj  =  ^ (arctg Rj+ arctg Rp  +  nj  ,  where  n^  =  -1,0;  for  ^€(-oo.^P  , 


n  =  0 


L(w)  =  -Ljj  -  iLj  =  U  -  4^1^^  U  U  -  cxpOjiYj  +  itry^  +  iJtYj  + 

Q(v 


i  arctg  R^)  , 


From  the  ctmdition  (3.1) 


Yj  *  i  (arctg  R^- arctg  rP +  113  ,  where  n^  =  -1 
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01  =  02  =  n3  =  0  from  the  condition  of  smoothness  of  the  solution  in  the  case  of  non- 
perforated  walls  (R  =  0). 

On  the  section  ^'e  [1,  d]  ,  t]  j  =  0  ,  d  =  exp  (•^)  in  the  plane  C  =  (see 

Fig  3b)  the  function  L(^(w  (z)))  in  the  general  case  has  a  complex  value 


UC(w(z)))  =  ILIcxpdtp^^)  .  (4) 

where  <Pl  =  9lo  +  <P.  • 

9^^  =  ARG[(iVf’  -  +  ARG[(i^/^)^2]  +  ARG[(iV^-  C/i]  . 


(b)  Let  us  find  an  analytical  function  F(^)  such,  that  in  the  plane 
^  =  exp(^  z)  =  -w^  =  ^'  +  iT|’  ,  having  two  cuts  along  the  real  axis  ^'e  (-«><»,  0)  and 
^'6  ^  1,  d  =  exp(-^)  J  ,  it  satisfies  the  following 

^'€  (-00, 0)  11'  =  ±0,  Im  F  =  0  , 

(5) 

^'6[l.d]  n’  =  ±0,ImF  = -9l  . 


In  general  form  the  function  F(Q  which  satisfies  condition  (5),  assumes  the  form: 


1  Va-k,  Va-k-  Va-k, 

IC  VC-1)  ^(d-0  ^ 


k.-’/a  lL,-'/a  lu-'/a 

^  (t-1)^  (d-t)^ 


+  P0+P1C 


where  ki  =  0, 1 ;  k2  =  0, 1 ;  k3  =  0, 1 ;  Pq  =  const ,  pj  =  const . 

From  this  it  will  be  seen  that  for  no  values  of  kj  can  the  condition  of  boundedness 
at  all  the  points  indicated  be  satisfied,  as  a  minimum  at  one  point  of  the  four  indicated 
is  not  limited,  unless  an  additional  condition  is  imposed  on  9l  for 
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Taking  this  into  account  we  demonstrate  that  with  the  condition  (6)  for  the  integer 
positive  values  kj  ,  k^,  k3  ,  such  that  k^  +  k^  +  k^  =  [  ® ,  Pp  =  =  0  ,  the  following 

equation  is  satisfied 


Vd^ 


a 


VTVd-t  Vt- 1  (t-Q 


Vs-k,  Vs-k,  V2-IU 

=  ;  ^;-i)  ^(d-C)  ^ 


a 

I  (Vl/  9.)  t' 


k,-'/2  k,-y2  k--'/:  dt 

^  (t-lp  (d-tp  —.(7) 

t-C 


It  is  sufficient  to  prove  the  correctness  of  equation  (7)  for  the  case  k]  =  k2  =  0  , 
k3  =  1 ,  since  the  remaining  ones  are  similar. 

Let  us  examine  the  difference  of  the  left-hand  and  right-hand  parts  in  expression  (7): 


which  in  fact  had  to  be  demonstrated.  Thus,  on  the  condition  (6),  the  required  analyticr^ 
function  assumes  the  form 
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(c)  Let  us  examine  the  analytical  function  T(0  =  L(Q  exp  (F(0)  .  The  function 

T(Q  has  discontinuities  along  the  line  0)  •  t)'  =  0  and  on  the  section 

^'€  [1,  d]  ,  Ti'  =  0  .  It  is  easy  to  see;  T(Q  =  Tj^  -  iTj  satisfies  conditions  on  the  line 
(-00, 0) .  ti'  =  0  which  are  analogous  to  the  boundary  condition  (3)  for  L(0  .  while 
on  the  section  ^'^[l.d],  Ti'  =  0  the  function  T(0  is  real. 

Let  us  introduce  a  further  subsidiary  function  <!>'  =  <I>/T ;  then  taking  into  account 
stages  (a)  and  (b)  for  d»’  the  boundary  conditions  in  the  plane  ^  assume  the  form: 


Ti’=±0,  Im«l>-  =  0;  1 

,  ^•€[l.d].  ImO)**  =  -h-(^)yT^(^)  ,  1 


ti’  =  ±0 


where  the  sign  «+»  corresponds  to  the  upper  edge  of  the  slit,  and  the  sign  «-» 
corresponds  to  the  lower  one. 

Let  us  represent  the  analytical  function  O'  as  the  sum  of  the  symmetrical  and  anti- 

symmetrical  analytical  functions:  O'  *  O'  +  o'  such,  that 

S  A 


OjC^’.tT)  =  Ug-iVj,  Vg(^',TT)  =  , 

I  I 

Thus  for  and  condition  (8)  will  give 

T]'  =  ±0 ,  ^'€  (-«,  0) ,  Im  O^*  =  0  ; 

I  Ti’  =  ±0,  ^'ell,d],  =  (^  +  ^)  ’/2  ; 
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I  T1’  =  ±0.  ^’€(^,0).  Im<I>3±  =  0; 

iT  =  ±0.  5'e[l.d].  Iin(D‘^±  =  (~  +  ~)  ’/2  . 

Taking  into  account  the  Zhukovskii  condition  on  the  trailing  edge  of  the  profile  and 
the  boundedness  of  the  function  to  infinity,  in  accordance  with  [7],  we  obtain 


d 


The  general  form  of  the  solution  for  O  assumes  the  form 


where  M„(C)  is  a  polynomia]  with  real  coefficients.  It  is  easy  to  sec  that  the  limit  solution 
(<I>  0 ,  X  ±09)  may  be  obtained  solely  on  condition  that  I  w)/Q(w)  =  (Ai  1)”^ 
and  Mn(0  =  Cq  =  const .  The  magnitude  of  the  constant  Cq  is  determined  from  the 
condition  of  absence  of  a  pole  at  the  point  ^  »  -Va  . 

The  associated  complex  velocity  may  be  writtm  in  the  following  femn: 
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<iKC(z))  =  u - iv  =  —  n V; - ^  (i Vo "  (i VC - 4.) ' 


X  exp 


a 

•  i-V^Vd-^  V^-  1  I  (<Pl  +  <P, 

71  J  o 


) 


dt 


Vt  Vd -  t  Vt-  1  (t -  0 


f  h'*^/exp(F^)  -  h  /exp(F  )  ^ 

J  lL(l)lcxp(i(<pL +‘P.))(t-C) 


>  . 


d 

,  I  d-C  f  {  h*/exp(F^)  -H  h7exp(F')^VrT _ ^ 

^  ^  ^  J  I  L(t)  I  cxp(i  (<Pl  -!  (p,))  Vd^t  (t  -  0 

where  C  =  "^’‘^(§*2)  ’ 


f(<p  +  q>,)  ds 


7t 


j  vivr=iv^(s-t) 


I 


=  ±Fo"‘(‘PLq'^‘P*^  ’  tell.d]  . 
for  the  following  conditions; 
d 


(9.1) 


f  h'*‘/cxp(-  ^q)  -  h'/expCPp) 

J  mt)l(i 


di  +  sign  (A)  X 


t  +  vl  A( 


-1/,  /  d  +  VlAf^  f y/exp(F-FQ)h7exp(FQ)) Vt^ 

•y  -/i  A  r*  + 1 1  I  L(t)  I  VTVd  - 1  (t + Vi  AT  rS 


dt  +  CQ  =  0,  (9.2) 
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h‘*’/exp(-FQ)  -  h  /cxp(Fq) 

IL(t)l(t  +  4]) 


dt  + 


h'^/cxp(-FQ)  +  h  /cxpOF^) 
IL(t)lVrVd^(i  +  ^j) 


Vi-  1  dt  +  Cq 


=  0  . 


(9.3) 


d 


J 


h'*'/cxp(-FQ)  -  h  /exp(FQ) 

I  L(t)  \(^) 


dt- 


+ 


h'*’/exp(-FQ)  +  h  /exp(F^ 
I  L(t)  I  VTVd-t  (t  +  ^2) 


Vt  -  1  dt 


+ 


=  0  . 


(9.4) 


Fig  4 
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Conditions  (9.3)  and  (9.4)  for  the  determination  of  Xj  and  X2  automatically  follow 

from  the  circumstance  that  during  passage  along  the  wall  through  the  points  Xj  and  X2  the 

pressure  coefficient  and  the  vertical  component  of  velocity  change  their  sign. 

In  essence,  solution  of  the  problem  is  reduced  to  determining  the  four  constants  A  , 
CO ,  xj ,  X2  .  A  diagram  of  the  solution  of  the  system  of  equations  (9.1)  to  (9.4)  is  given 
by  Fig  4. 

The  pressure  coefficient  on  the  porous  boundary  is  equal  to 


r  d 


h’^/expC-r^o)  -  h~/exp(Fo) 
mt)l(t-0 


dt± 


d 


h'^/exp(-Fo)  -  h'/exp(Fp) 
I  L(t)  I  VrVd^(t- 0 


Vt-  1  dt Cp  > 


where 


cos  (-arctgR^) 

cos  (-^uctg  Rt  +  rry.) 

G=  <  ,  . 

cos  (-arctg  +  7t(Ti  +  72^ ' 

cos  (-arctg  R^  +  7t(7j  +72  +  73)) 
C=  Lowcrwall 

C=  Upper  wall 

\  2/x>x, 


The  pressure  coefficient  on  the  profile  is  equal  to 
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c*  =  -2U*  =  i 

P  71 


L(C(x  +  iO))  X 


/  “ 

X  exp”^  tvcvdrcvm  J 


a 

I*  h’^/expC-pQ)  -  h'/expCF^) 

J  II 


((p  +  (p,)  dt 
^0 


VTVd-t  Vt- 1  (t-0 


X 


iL(t)i  vrvd-t(t-o 


/exp(-F  )  -  h  /expOF  )  ^ _ 

- - 2-V7rTdt  +  c  , 


where  the  sign  «+»  denotes  the  upper  surface  of  the  profile,  and  «-»  the  lowf  r. 


3  Calculations  of  the  pressure  on  the  profile,  which  were  performed  for  comparison  in 
the  simplest  case  (symmetrical  problem  of  boundary  conditions  on  the  upper  and  lower 
walls,  or  H  -4  «»),  gave  a  good  agreement  with  the  results  of  the  calculation  made  in 
Refs  1, 2  and  S. 

Figs  5  and  6  show  calculations  of  flow  about  a  plate  at  an  angle  of  attack  of  3°  in  a 
tunnel  (H  =  6)  dependent  upon  boundary  conditions  in  porous  walls. 

As  follows  from  Refs  3  and  4,  symmetrical  assigning  on  the  upper  and  lower  walls 
of  the  tunnel  of  the  parameter  of  porosity  R  corresponds  to  a  differing  coefficient  of 
porosity  of  these  walls. 

As  a  rule,  in  wind  tunnels,  the  coefficient  of  porosity  is  identical  for  the  upper  and 
lower  wall.  As  will  be  seen  from  the  calculations  performed,  tlie  non-symmetrical  assigning 
of  the  boundary  conditions  on  the  porous  surface  perceptibly  changes  the  flow  picture.  A 
change  over  the  length  of  the  tunnel  of  the  parameter  of  porosity  on  the  upper  wall  has  a 
slight  effect  on  the  distribution  of  pressure  on  the  boundary,  while  it  is  practically  absent  in 
the  region  of  location  of  the  plate. 

Therefore,  m  the  case  of  an  identical  coefficient  of  porosity  of  the  walls  it  is 
necessary  to  take  into  account  die  difference  in  the  parameter  of  porosity  dependent  upon  the 
direction  of  flow  of  the  gas,  that  is,  into  the  tunnel  or  out  of  it 

The  author  wishes  to  thank  V,M.  Neiland  and  O.K.  Semenov  for  useful 
discussions. 
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